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Abstract. We study the operator product expansion of the plaquette (gluon condensate) and the self-energy of an infinitely 
heavy quark. We first compute their perturbative expansions to order and respectively, in the lattice scheme. In both 
cases we reach the asymptotic regime where the renormalon behavior sets in. Subtracting the perturbative series, we obtain 
the leading non-perturbative corrections of their respective operator product expansions. In the first case we obtain the gluon 
condensate and in the second the binding energy of the heavy quark in the infinite mass limit. The results are fully consistent 
with the expectations from renormalons and the operator product expansion. 
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INTRODUCTION 

The operator product expansion (OPE) [1] is a fundamental tool for theoretical analyses in quantum field theories. Its 
validity is only proven rigorously within perturbation theory, to arbitrary finite orders [2]. The use of the OPE in a 
non-perturbative framework was initiated by the ITEP group [3] (see also the discussion in Ref. [4]), who postulated 
that the OPE of a correlator could be approximated by the following series: 

correlator)^) ~ ^ ^Cd{a){Od), (1) 

where the expectation values of local operators Oj are suppressed by inverse powers of a large external momentum 
Q ^ Aqcd. according to their dimensionality d. The Wilson coefficients Q(o:) encode the physics at momentum 
scales larger than Q. These are well approximated by perturbative expansions in the strong coupling parameter a: 

Q(a)~ ^c„a"+'. (2) 

H>0 

The large-distance physics is described by the matrix elements (Od) that usually have to be determined non- 
perturbatively: {Od) 

It can hardly be overemphasized that (except for direct predictions of non-perturbative lattice simulations, e.g., on 
light hadron masses) all QCD predictions are based on factorizations that are generalizations of the above generic 
OPE. 

There exist some major questions related to the OPE that have to be addressed:^ 

• Are the perturbative expansions of Wilson coefficients asymptotic series? 


* There are also some important points that we do not address here, including: 

We do not consider ambiguities associated to short-distance non-perturbative effects, which would give rise to singularities further away from the 
origin of the Borel plane than those we study here. 

We take the validity of the OPE in pure perturbation theory for granted. This assumption is solid in cases with a single large scale, Q^, and in 
Euclidean spacetime. 

We will not discuss the validity of the (non-perturbative) OPE for timelike distances that can occur in Minkowski spacetime, an issue related to 
possible violations of quark-hadron duality. 



• If so: are the associated ambiguities of the asymptotic behavior consistent with the OPE, i.e. with the positions of 
the expected renormalons [5] in the Borel plane? 

• Is the OPE valid beyond perturbation theory? 

• What is the real size of the first non-perturbative correction within a given OPE expansion? 

• Is this value strongly affected by ambiguities associated to renormalons? 

In this paper we summarize and discuss our recent results [6-10], which address these questions for the case of 
the plaquette and the energy of an infinitely heavy quark in the pure gluodynamics approximation to QCD. Both 
analyses utilize lattice regularization. Contrary to, e.g., dimensional regulation, lattice regularization can be defined 
non-perturbatively. Using a lattice scheme rather than the MS scheme, we can, not only expand observables in 
perturbation theory, but also evaluate them non-perturbatively. Another advantage of this choice is that it enables us to 
use numerical stochastic perturbation theory [11-13] to obtain perturbative expansion coefficients. This allows us to 
realize much higher orders than would have been possible with diagrammatic techniques. A disadvantage of the lattice 
scheme is that, at least in our discretization, lattice perturbative expansions converge slower than expansions in the MS 
coupling. This means that we have to go to comparatively higher orders to become sensitive to the asymptotic behavior. 
Many of the results obtained in a lattice scheme either directly apply to the MS scheme too or can subsequently easily 
(and in some cases exactly) be converted into this scheme. 

In our studies we used the Wilson gauge action [14]. We define the vacuum expectation value of a generic operator 
B of engineering dimension zero as 

{B) = (a|B|f2) = (3) 

with the partition function Z = and measure [dUx,fi] = 1^) denotes the vacuum state, 

Ae is a Euclidean spacetime lattice with lattice spacing a and Ux.fi ~ g SU(3) is a gauge link. 


THE PLAQUETTE: OPE IN PERTURBATION THEORY 


Eor the case of the plaquette we have B P, where 
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and Ux.fxv denotes the oriented product of gauge links enclosing an elementary square (plaquette) in the /r-V plane of 
the lattice. Eor details on the notation and simulation set-up see Ref. [9]. 

(P) will depend on the lattice extent Na, the spacing a and a = / (47r) = a(a^*) (note that a is the bare lattice 

coupling and its natural scale is of order a^^). We first compute this expectation value in strict perturbation theory. 
In other words, we Taylor expand in powers of g before averaging over the gauge configurations (which we do using 
NSPT [11-13]). The outcome is a power series in a: 




NSPT 


n>0 


The dimensionless coefficients Pn{N) are functions of the linear lattice size N. We emphasize that they do not depend 
on the lattice spacing a or on the physical lattice extent Na alone but only on the ratio N = {Na)/a. 

We are interested in the large-A (i.e. infinite volume) limit. In this situation 


1 1 
a Na 


(5) 


and it makes sense to factorize the contributions of the different scales within the OPE framework. The hard modes, 
of scale ~ 1/a, determine the Wilson coefficients, whereas the soft modes, of scale ~ 1/ {Na), can be described by ex¬ 
pectation values of local gauge invariant operators. There are no such operators of dimension two. The renormalization 
group invariant definition of the gluon condensate 
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is the only local gauge invariant expectation value of an operator of dimension in pure gluodynamics. In the purely 

perturbative case discussed here, this only depends on the soft scale 1 / (Na), i.e. on the lattice extent. On dimensional 
grounds, the perturbative gluon condensate (G^)soft is proportional to 1/ {Na)^, and the logarithmic (A^a)-dependence 
is encoded in a[l/{Na)]. Therefore, 

= -^ E /„a"+'[l/(iVa)], ( 7 ) 


and the perturbative expansion of the plaquette on a finite volume of sites can be written as 


( 7 ’)pert(iV) — 7 pert( 0 !)(l) + 2 ^CG(o:)a (G + ) ( 8 ) 

where 

7pert(a) = (^9) 

n>0 

and pn are the infinite volume coefficients that we are interested in. The constant pre-factor n'^/36 is chosen such that 
the Wilson coefficient, which only depends on a, is normalized to unity for a = 0. It can be expanded in a: 

CG(a) = l + Ec-t«'+'- (10) 

k>0 


Since our action is proportional to the plaquette P, Cq is fixed by the conformal trace anomaly [15, 16]: 
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The j3-function coefficients^ Pj are known in the lattice scheme for j <3 (see Eq. (25) of Ref. [9]). 
Combining Eqs. (7), (8) and (10) gives 


(R)pert(A^) — E 
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where fn{N) is a polynomial in powers of In(A^). Eitting this equation to the perturbative lattice results, the first 35 
coefficients were determined in Ref. [9]. The results were confronted with the expectations from renormalons: 
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where b = pi/{2p§), si = {pf - j3oj32)/(4hj3^) and S 2 = {Pf - 2 .popip 2 + PqP 3 )/{^Bb^p^) are defined so that 
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In Eig. 1 we compare the infinite volume ratios p„/(npn-i) to the expectation Eq. (14): the asymptotic behavior of the 
perturbative series due to renormalons is reached around orders n ~ 27 — 30, proving, for the first time, the existence 
of the renormalon in the plaquette. Note that incorporating finite volume effects is compulsory to see this behavior, 
since there are no infrared renormalons on a finite lattice. To parameterize finite size effects we made use of the purely 
perturbative OPE Eq. (8). The behavior seen in Eig. 1, although computed from perturbative expansion coefficients, 
goes beyond the purely perturbative OPE since it predicts the position of a non-perturbative object in the Borel plane. 


^ We define the /1-function as ji (a) = da j din p = —/lo/(2;r)a^ ^ Pi/ (8n:^)a^ — ■ ■ ■, ie. j3o = 11. 
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FIGURE 1. The ratios pn/{npn-i) compared with the leading order (LO), next-to-leading order (NLO), NNLO and NNNLO 
predictions of the 1/n-expansion Eq. (14). Only the ‘W = =” extrapolation includes the systematic uncertainties. We also show 
finite volume data for N = 28, and the result from the alternative N ^ extrapolation including some l/N^ corrections. The 
symbols have been shifted slightly horizontally. 


THE PLAQUETTE: OPE BEYOND PERTURBATION THEORY 


Since in NSPT we Taylor expand in powers of g before averaging over the gauge variables, no mass gap is generated. 
In non-perturbative Monte-Carlo (MC) lattice simulations an additional scale, Aqcd ^ is generated 

dynamically (see also Eq. (15)). However, we can always tune N and a such that 
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In this small-volume situation we encounter a double expansion in powers of a/{No) and aAqcd [or, equivalently, 
(Aa)AQCD X a/{No)}. The construction of the OPE is completely analogous to that of the previous section and we 
obtain^ 


{P)mc = ^ 


MC 


— Ppert(ct)(l) + ^C’G(o:)a (G )mc + ^(^ )■ 


(17) 


In the last equality we have factored out the hard scale 1/a from the scales 1 /{Na) and Aqcd, which are encoded in 
(G2) MC- Exploiting the right-most inequality of Eq. (16), we can expand (G^)mc as follows: 


(g2)mc = {1 + ^[A2cD(Aa)2]} . 


(18) 


^ In the last equality, we approximate the Wilson coefficients by their perturbative expansions, neglecting the possibility of non-perturbative 
contributions associated to the hard scale l/a. These would be suppressed by factors ~ exp(—2;r/a) and therefore would be sub-leading, relative 
to the gluon condensate. 















Hence, a non-perturbative small-volume simulation would yield the same expression as NSPT, up to non-perturbative 
corrections that can be made arbitrarily small by reducing a and therefore Na, keeping N fixed. In other words, 
^NSPT(^) _ jQ non-perturbative corrections. 

We can also consider the limit 

-»Aqcd>—• (19) 

a Na 

This is the standard situation realized in non-perturbative lattice simulations. Again the OPE can be constructed as in 
the previous section, Eq. (17) holds, and the and c„-values are still the same. The difference is that now 

(G^)mc = (G^)np 1 + ^ ( — . '2 J , (20) 

[ V'^QCd(^«) / J 

where (G^)np ~ ^qcd so-called non-perturbative gluon condensate introduced in Ref. [3]. Erom now on we will 
call this quantity simply the “gluon condensate” (G^). We are now in the position 

• to determine the gluon condensate and 

• to check the validity of the OPE (at low orders in the scale expansion) for the case of the plaquette. 

In order to do so we proceed as follows. The perturbative series is divergent due to renormalons and other, sub-leading, 
instabilities."^ This makes any determination of (G^) ambiguous, unless we define precisely how to truncate or how to 
approximate the perturbative series. A reasonable definition that is consistent with (G^) ~ Aq^p can only be given if 
the asymptotic behavior of the perturbative series is under control. This has only been achieved recently [9], where the 
perturbative expansion of the plaquette was computed up to see the previous section. The observed asymptotic 

behavior was in full compliance with renormalon expectations, with successive contributions starting to diverge for 
orders around within the range of couplings a typically employed in present-day lattice simulations. 

Extracting the gluon condensate from the average plaquette was pioneered in Refs. [17-20] and many attempts 
followed during the next decades, see, e.g.. Refs. [21-30]. These suffered from insufficiently high perturbative orders 
and, in some cases, also finite volume effects. The failure to make contact to the asymptotic regime prevented a reliable 
lattice determination of (G^). This problem was solved in Ref. [10], which we now summarize. 

Truncating the infinite sum at the order of the minimal contribution provides one definition of the perturbative series. 
Varying the truncation order will result in changes of size Aq^^qa"^, where the dimension c/ = 4 is fixed by that of the 
gluon condensate. We approximate the asymptotic series by the truncated sum 

n 

Sp{a) = S„g{a), where 5„(a) = ^ (21) 

./=o 

no = no(a) is the order for which is minimal. We then obtain the gluon condensate from the relation 

(G^) = [(P)Mc(a) -5Ka)] + . (22) 

Cq* (a) is proportional to the j3-function, and the first few terms are known, see Eq. (11). The corrections to Cg = 1 
are small. However, the and ^(a^) terms are of similar sizes. We will account for this uncertainty in our error 

budget. 

Eollowing Eq. (22), we subtract the truncated sum Sp{a) calculated from the coefficients pn of Ref. [9] from the MC 
data on {P)mc{oc) of Ref. [31] in the range j3 G [5.8,6.65] (j3 = where fl(j3) is given by the phenomenological 
parametrization of Ref. [32] (x = j3 — 6) 

a = ro exp (-1.6804 - 1.733 lx -f 0.7849x^ - 0.4428x^) , (23) 

where ro « 0.5 fm. This corresponds to (a/ro)"^ G [3.1 x 10^^,5.5x 10^^], covering more than two orders of magnitude. 


The leading renormalon is located atM = ^/2 = 2in the Borel plane, while the first instanton-anti-instanton contribution occurs at w = j3o = 

1Wg/3-11»2. 





FIGURE 2. Left panel: Eq. (22) evaluated using the = 16 and N = 32 MC data of Ref. [31]. The N = 32 outer error bars 
include the error of Sp(a). The error band is our prediction for {Oq), Eq. (24). 

Right panel: Eq. (22) times o' vs. fl^(a)/ry from Eq. (23). The linear fit with slope Eq. (24) is to the < 0.0013 Tq points only. 


Multiplying this difference by (TT^Coa^) gives rQ{OQ) plus higher order non-perturbative terms. We show this 
combination in the left panel of Fig. 2. The smaller error bars represent the errors of the MC data, the outer error bars 
(not plotted for = 16) the total uncertainty, including that of Sp. This part of the error is correlated between different 
j3-values (see the discussion in Ref [10]). The MC data were obtained on volumes = 16^ and = 2>2‘^. Towards 
large p -values the physical volumes [A^a(j3 will become small, resulting in transitions into the deconhned phase. For 
P < 6.3 we hnd no significant differences between the = 16 and N = 32 results. In the analysis we restrict ourselves 
to the more precise N = 32 data and, to keep finite size effects under control, to j3 < 6.65. We also limit ourselves 
to P > 5.8 to avoid large GicP') corrections. At very large )3-values not only does the parametrization Eq. (23) break 
down but obtaining meaningful results becomes challenging numerically: the individual errors both of (F’)mc(oc) 
of Sp{a) somewhat decrease with increasing p. However, there are strong cancellations between these two terms, in 
particular at large j3-values, since this difference decreases with exp{l6n^P/33) on dimensional grounds 

while (P)mc depends only logarithmically on a. 

The data in the left panel of Fig. 2 show an approximately constant behavior.^ This indicates that, after subtracting 
Sp{a) from the corresponding MC values (F)mc(o!), the remainder scales like a^. This can be seen more explicitly 
in the right panel of Fig. 2, where we plot this difference in lattice units against a^. The result is consistent with a 
linear behavior but a small curvature seems to be present that can be parametrized as an a®-correction. The right-most 
point (P = 5.8) corresponds to ~ 1.45 GeV while p = 6.65 corresponds to ~ 5.3 GeV. Note that a^-terms are 
clearly ruled out. 

We now determine the gluon condensate. We obtain the central value and its statistical error (G^) = 3.177(36)rQ ^ 
from averaging the N = 32 data for 6.0 <P<6 .65. We now estimate the systematic uncertainties. Different infinite 
volume extrapolations of the Pn{N) data [9] result in changes of the prediction of about 6%. Another 6% error is due to 
including an a^-term or not and varying the fit range. Next there is a scale error of about 2.5%, translating into units 
of KQ. The uncertainty of the perturbatively determined Wilson coefficient Cq is of a similar size. This is estimated 
as the difference between evaluating Eq. (11) to 0'{a?') and to Adding all these sources of uncertainty in 

quadrature and using the pure gluodynamics value [33] A^ = 0.602(48)rQ^^ yields 

(G^) = 3.18(29)ro 4 = 24.2(8.0)A^. (24) 

The gluon condensate Eq. (6) is independent of the renormalization scale. However, (G^) was obtained employing 
one particular prescription in terms of the observable and our choice of how to truncate the perturbative series within a 
given renormalization scheme. Different (reasonable) prescriptions can in principle give different results. One may for 


^ Note that increases from 26 to 27 at j3 = 5.85, from 27 to 28 at j3 = 6.1 and from 28 to 29 at j3 = 6.55. This quantization of no explains the 
visible jump at j3 =6.1. 
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FIGURE 3. Differences {P)mc{^) — Sn{(x) between MC data and sums truncated at orders a"+* (5_i = 0) vs. a{a)/rQ. The 
lines cx: flt are drawn to guide the eye. 


instance choose to truncate the sum at orders no{a) ± yno(^ and the result would still scale like We estimated 

this intrinsic ambiguity of the definition of the gluon condensate in Ref. [9] as 5{G^) = 36/{n^Ccci'^) 
i.e. as ^no(a) times the contribution of the minimal term: 

5(g2)=27(11)A|^. (25) 

Up to 1 /no-corrections this definition is scheme- and scale-independent and corresponds to the (ambiguous) imaginary 
part of the Borel integral times 

In QCD with sea quarks the OPE of the average plaquette or of the Adler function will receive additional contribu¬ 
tions from the chiral condensate. For instance (G^) needs to be redefined, adding terms (7„,(a)my/y/) [34, 35]. Due 
to this and the problem of setting a physical scale in pure gluodynamics, it is difficult to assess the precise numerical 
impact of including sea quarks onto our estimates 

(G^) ~ 0.077 GeV'^, 5 (G^) ~ 0.087 GeV'^, (26) 

which we obtain using ro ~ 0.5 fm [36]. While the systematics of applying Eqs. (24)-(25) to full QCD are unknown, 
our main observations should still extend to this case. We remark that our prediction of the gluon condensate Eq. (26) 
is significantly bigger than values obtained in one- and two-loop sum rule analyses, ranging from 0.01 GeV^ [3, 37] 
up to 0.02 GeV^ [38, 39]. However, these numbers were not extracted in the asymptotic regime, which for a d = A 
renormalon in the MS scheme we expect to set in at orders n>l. Moreover, we remark that in schemes without a hard 
ultraviolet cut-off, like dimensional regularization, the extraction of (G^) can become obscured by the possibility of 
ultraviolet renormalons. Independent of these considerations, all these values are smaller than the intrinsic prescription 
dependence Eq. (25). 

Our analysis confirms the validity of the OPE beyond perturbation theory for the case of the plaquette. Our de¬ 
scaling clearly disfavors suggestions about the existence of dimension two condensates beyond the standard OPE 
framework [25, 40-43]. In fact we can also explain why an a^-contribution to the plaquette was found in Ref. [25]. 
In the log-log plot of Fig. 3 we subtract sums truncated at different fixed orders from (P)mc- The scaling 

continuously moves from ^ at to ^ d around Note that truncating at an a-independent fixed order 














is inconsistent, explaining why in the figure we never exactly obtain an a^-slope. For n ~ 9 we reproduce the a^-scaling 
reported in Ref. [25] for a fixed order truncation at n = 7. In view of Fig. 3, we conclude that the observation of this 
scaling power was accidental. 


THE BINDING ENERGY OF HQET: OPE IN PERTURBATION THEORY 

The OPE of the plaquette is analogous to the OPE of the vacuum polarization (or the Adler function) in position 
space. As we already mentioned, the OPE concept of the factorization of scales can also be applied to more general 
kinematical settings (in particular to cases where some scales are defined in Minkowski spacetime). A prominent 
example is heavy quark effective theory (HQET). In this case the term (Oi) of Eq. (1) is replaced by a non-perturbative 
quantity, the so-called heavy quark binding energy A, that cannot be represented as an expectation value of a local 
gauge invariant operator. We consider the self-energy of a heavy quark in the infinite mass limit (in other words, the 
self-energy of a static source in the triplet representation, for other color representations see Ref. [7]), that is closely 
related to A. We compute this quantity in close analogy to the case of the plaquette in lattice regularization. 

Eirst we compute the self-energy of the static quark in perturbation theory (again using NSPT). We obtain this from 
the Polyakov loop in an asymmetric volume N^Njo^ of spatial and temporal extents Nsa and Nto, respectively: 

1 1 Ni~l 

L{Ns,Nt) = ^Y.^T^^ n , (27) 

X ^ U4/a=0 

or, more specifically, from its logarithm 

PiNs,Nr) = ^ . (28) 

aJyr „>o 

Again, denotes a gauge link and x = (x,X 4 ) G Ae are Euclidean lattice points. We define the energy of a static 
source and its perturbative expansion in a finite spatial volume, 

dm{Ns)= lim P(A(s,Ar) = - Y c„(As)a"+'(l/a), where c„{Ns) = lim c„{Ns,Nt) , (29) 

Nj-^oo a Nj^oo 


and its infinite volume limit 

1 °° 

5m = lim P(A 5 ,Ar) = - Y c„a"^Vl/a) with c„ = limc„(A( 5 ). ( 

We now construct the purely perturbative OPE in a finite volume. Eor large Ns, we can write [a^* ^ 

4 1 fn{Ns)\ „,n+Pi ix , 1 \ . 


5m{Ns) = - Y 


a"+'(l/fl)-l-^ 


««>o 


Note the similarity between this equation and Eq. (12). fn{Ns) = is again a polynomial in powers 

of In (As), and the ^ are known combinations of j3-function coefficients and lower order infinite volume expansion 
coefficients Ck, k <n. The main difference with respect to the gluon condensate is that the power correction scales like 
1 /As, instead of 1 /A^, and that now the Wilson coefficient is trivial. This 1 /As scaling also means that the renormalon 
behavior will show up at lower orders n of the perturbative expansion. Eitting the Cn{Ns) data to this equation, the first 
20 Cn coefficients were determined in Refs. [6-8] and confronted with the renormalon expectations:® 

V r(«-|-1 +b) , ^.^1 b^{s\/ 2 -S'^ _ 

" \ 2 ;r/ r(l-l-fi) {n + b) {n + b){n + b—l) 


Here we deviate from Refs. [6-8] in the definition of the constant S 2 , see Eq. (15). 



9 
8 
7 

^ 6 

£Z 

S 5 
4 

3 
2 
1 

FIGURE 4. Ratios c„/(nc„_i) of the unsmeared (squares) and smeared (circles) triplet static self-energy coefficients c„ in 
comparison to the theoretical prediction Eq. (33), truncated at different orders in 1/n. In the inset we magnify the asymptotic 
region. 


smeared 

LO 


I—•—I 



J_L 


8 10 12 14 16 18 


J_I_L 


2 4 6 8 10 12 14 16 18 

n 


In the lattice scheme the numerical values of the above coefficients read bs\ = 1.36095381(11) and b^{s\/2 — S 2 ) = 
5.34(51). As expected, the above expansion converges much faster in 1/n than Eq. (13). Calculating the ratio of 
subsequent perturbative coefficients gives 
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(33) 


We remark that Eq. (14) includes the effect of the non-trivial Wilson coefficient Cq. Therefore, just setting r/ = 1 in 
that equation does not result in Eq. (33) above. 

In Eig. 4 we compare the data to Eq. (33) for two different lattice discretizations of the covariant temporal derivative, 
which amounts to “smearing” or not smearing temporal gauge links. This should not affect the infrared behavior and, 
indeed, beyond the first few orders the difference becomes invisible. The asymptotic behavior of the perturbative series 
due to the renormalon is confirmed in full glory for n > 8, proving the existence of the renormalon behavior in QCD 
beyond any reasonable doubt. Again the incorporation of finite volume effects was decisive to obtain this result. 


THE BINDING ENERGY OF HQET: OPE BEYOND PERTURBATION THEORY 

The methods used for the gluon condensate can also be applied to other observables. We now consider the non- 
perturbative evaluation of the energy of a static-light meson on the lattice and its OPE: 

£^Mc(0 !) = £^pert(0!) + A+^(aAQ(;;p) . (34) 

A is the non-perturbative binding energy and £’pert(oc) = 5m{a) is the self-energy of the static source in perturbation 
theory, i.e. Eq. (30). In HQET the mass of the B meson is given as ms = '«fo,pert + A+ ^(l/'w*), where m^ pert is the 
l7-quarkpole mass which suffers from the same renormalon ambiguity as 5m [44^6]. 
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FIGURE 5. aA = oEmc — o^Sgm vs. fl/ro- The sum aSgm Eq. (36) was also converted into the MS scheme at two (MS 2 ) and three 
(MS 3 ) loops and truncated at the respective minimal orders. The curves are fits to Aa + ca^ frQ. 



The perturbative expansion of adm{a) — was obtained in Refs. [ 6 , 7] up to see the previous 

section. Its intrinsic ambiguity 

5A = v^c„oa"o+' = 0.748(42 )Ams = 0.450(44)ro ^ (35) 

was computed in Refs. [7, 8 ]. MC data for the ground state energy Emc of a static-light meson with the Wilson action 
can be found in Refs. [47^9]. 

Truncating the infinite sum in 5m at the order of the minimal contribution provides one definition of the perturbative 
series. As we did for the plaquette, we approximate the asymptotic series by the truncated sum 

1 "o(«) 

Ss„,{a) = -Y. > (36) 

« j=0 

where no = no{a) is the order for which c„Qa"o+* is minimal. While for the gluon condensate we expected an a^- 
scaling (see the right panel of Fig. 2), for aA = uEmc — ciSsm{(^) we expect a scaling linear in a. Comforting enough 
this is what we find, up to the expected aff{aAQ^^) discretization corrections, see Fig. 5. Subtracting the partial sum 

truncated at orders no (a) =6 from the j3 S [5.9,6.4] data, we obtain A = 1.55(8)rQ ^ from such a linear plus quadratic 
fit, where we only give the statistical uncertainty. The errors of the perturbative coefficients are all tiny, which allows 
us to transform the expansion a5m{a) into MS-like schemes and to compute A accordingly. We define the schemes 
MS 2 and MS 3 by truncating = a(l + d]_a + d 20 p' -I-...) exactly at and respectively. The 

dj are known for j <3 [7, 8 ]. We typically find n^®'(aj^ ) =2,3 and obtain A ^ 2.17(8)rQ * and A ^ 1.89(8)rQ ^ 
respectively, see Fig. 5. We conclude that the changes due to these resummations are indeed of the size 5A ^ 0.5rg *, 
adding confidence that our definition of the ambiguity Eq. (35) is neither a gross overestimate nor an underestimate. 
For the plaquette, where we expect ^ 7, we cannot carry out a similar analysis, due to the extremely high precision 
that is required to resolve the differences between Sp{a) and (7 ’)mc(<t), which largely cancel in Eq. (22). 






CONCLUSIONS 


For the first time ever, perturbative expansions at orders where the asymptotic regime is reached have been obtained 
and subtracted from non-perturbative Monte Carlo data of the static-light meson mass and of the plaquette, thereby 
validating the OPE for these cases beyond perturbation theory. The scaling of the latter difference with the lattice 
spacing confirms the dimension d = A. Dimension d < A slopes appear only when subtracting the perturbative series 
truncated at fixed pre-asymptotic orders. Therefore, we interpret the lower dimensional “condensates” discussed in 
Ref. [40] as approximate parametrizations of unaccounted perturbative effects, i.e. of the short-distance behavior. 
These will be observable-dependent, unlike the non-perturbative gluon condensate. Such simplified parametrizations 
introduce unquantifiable errors and, therefore, are of limited phenomenological use. As we have demonstrated above 
(see Fig. 3), even the effective dimension of such a “condensate” varies when truncating a perturbative series at 
different orders. In Refs. [41^3] various analyses, based on models such as string/gauge duality or Regge models, 
have been made claiming the existence of non-perturbative dimension two condensates. Our results strongly suggest 
that there may be flaws in these derivations. 

We observe that A and (G^) do not depend on the lattice spacing (i.e. on the renormalization scale). In other words, 
they are renormalization group invariant quantities, as expected. This is coherent with the interpretation of these 
quantities to be of order Aqcd and Aq^^, respectively. However, the values of A and (G^) will depend on the details 
of how the divergent perturbative series is truncated or estimated (and therefore implicitly also on the scheme used in 
the perturbative expansion), as well as on the observable used as an input in the determination. 

We have obtained an accurate value of the gluon condensate in SU(3) gluodynamics, Eq. (24). It is of a similar 
size as the intrinsic difference, Eq. (25), between (reasonable) subtraction prescriptions. This result contradicts the 
implicit assumption of sum rule analyses that the renormalon ambiguity is much smaller than leading non-perturbative 
corrections. The value of the gluon condensate obtained with sum rules can vary significantly due to this intrinsic 
ambiguity if determined using different prescriptions or truncating at different orders in perturbation theory. Clearly, 
the impact of this, e.g., on determinations of a, from T-decays or from lattice simulations needs to be assessed 
carefully. 

As already mentioned in the previous paragraph, due to the non-convergent nature of the perturbative series, the 
binding energy A and the gluon condensate (G^) are in principle ill-defined quantities. The intrinsic ambiguity of 
their definition can, however, be estimated. The ambiguity of the HQET binding energy 5A = 0.75 (4) A-^ as well as 
the ambiguity 5{G^) = 27(11)A^ of the non-perturbative gluon condensate are scale- and renormalization-scheme 
independent, at least up to 1 /no-corrections, where no is the order of the minimal term of the perturbative series. 
In the first case this ambiguity is significantly smaller than A ~ 2.6A^. In the second case the ambiguity is larger 
than values typically quoted for (G^), including our result (G^) ~ 24A^. The size of 5{G^) > (G^) means that, 
at least in pure gluodynamics,^ (G^) should not be used to estimate the magnitude of unknown non-perturbative 
corrections. Instead, the ambiguity 5{G^) should be used for this purpose. An exception to this rule are situations 
where the renormalon ambiguities cancel exactly or are small. For instance, the perturbative expansion of a difference 
between two observables {B 2 ) — {Bi) that receive contributions 0 = Cg (G^) jQ^ with a relative normalization such that 
Cbi = Cb 2 + ^{cc) will be partially blind to the associated d = A infrared renormalon if the (Bj) are expanded at the 
same scale, in the same renormalization scheme and the same method is used to truncate both perturbative series. 
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